Roles of dark energy perturbations in the dynamical dark energy models: 

Can we ignore them? 
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We show the importance of properly including the perturbations of the dark energy component 
in the dynamical dark energy models based on a scalar field and modified gravity theories in order 
to meet with present and future observational precisions. Based on a simple scaling scalar field 
dark energy model, we show that observationally distinguishable substantial differences appear by 
ignoring the dark energy perturbation. By ignoring it the perturbed system of equations becomes 
inconsistent and deviations in (gauge-invariant) power spectra depend on the gauge choice. 
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The high-z type la supernovae (SNIa) luminosity- 
distance relation suggests that the expansion rate of our 
universe is currently under acceleration [11]. The cos- 
mological constant is readily (re)introduced to explain 
the observation theoretically. Theoretical studies of the 
large-scale structure formation process imprinted in the 
matter power spectrum [2|, Q and the cosmic microwave 
background radiation (CMB) power spectrum Q also 
favor presence of substantial amount of agent with re- 
pulsive nature like the cosmological constant. With the 
advent of the recent acceleration, the long lasting age- 
problem of the world model, which has persisted ever 
since the first observation of the expansion of the uni- 
verse, has now evaporated from the cosmological scene. 

The nature of the agent causing the acceleration, how- 
ever, is still unknown and it is one of the fundamental 
mysteries in the present day theoretical cosmology. Al- 
though the cosmological constant is a historically well 
known possibility, it also has two well appreciated prob- 
lems: the cosmological constant (why so small) problem 
and the coincidence (why now or fine tuning) problem. 
Large amount of literature has been devoted to address 
these problems, especially the latter one, by introduc- 
ing dynamical agents, often termed the dark energy. As 
far as we can tell the fine tuning problem has not been 
properly addressed even using the dynamical dark en- 
ergy. Introduction of the dynamic possibility of the dark 
energy, however, has opened a whole new arena for cos- 
mological research based on variety of possibilities using 
field, fluid, modified gravity, other dimensions, etc. 

In the case of the cosmological constant as the dark en- 
ergy, due to its constant nature (both in time and space) 
its contribution directly appears only in the background 
world model. However, when we consider the dynami- 
cal dark energy we should pay attention to its dynamical 
roles not only in the background world model but also 
in the structure formation process. Here, we address the 
importance of properly including the role of dark energy 
perturbation (DEP) imprinted in the large-scale matter 
and the CMB anisotropics power spectra, and the pertur- 



bation growth process especially in the context of present 
and future observations with due precision. 

Recent dramatic progresses made in the observational 
cosmology open possibility to constrain the character of 
the dark energy, and call for equally precise theoretical 
tools in the cosmic structure formation process. The ex- 
pansion history based on the SNIa, the matter power 
spectrum, the CMB anisotropy power spectra, and the 
perturbation growth factor provide four domains where 
theories meet with observations. The relevant present 
and future observation programs in the CMB, SNIa, and 
the large-scale clustering include the WMAP (Wilkin- 
son Microwave Anisotropy Probe) and the Planck mis- 
sions, 2dFGRS (Two-degree-Field Galaxy Redshift Sur- 
vey), SDSS (Sloan Digital Sky Survey), to mention a few. 
The large-scale clustering can be probed by diverse obser- 
vations: weak lensing, Lyman-a and hydrogen 21cm to- 
mography. X-ray galaxy clustering mass function, galaxy 
redshift-space distortion, integrated Sachs- Wolfe effect, 
etc. In the following we will compare our results with 
SDSS DR7 (seventh Data Release) for the matter power 
spectrum g, WMAP 5-year data for the CMB spectrum 
[j| , and the future X-ray and weak lensing observations of 
clusters using X-ray surveys for the perturbation growth 
factor g. 

Our study is motivated by often used practices in the 
literature which ignore the DEP even in the case of dark 
energy models using the scalar field or modified gravity 
theories, see That is, in the presence of a dynami- 
cal dark energy it is not guaranteed to use the following 
conventionally known equation Q 



2HSb - AiiGSpb = 0, 



(1) 



which is true only for the cosmological constant as the 
dark energy; Sf, = Spb/pb is the relative density fluctu- 
ation of baryon component, H = a/a, a is the cosmic 
scale factor, and an overdot denotes a time derivative. 
In the presence of dynamical dark energy we have contri- 
butions from the DEP in the right-hand side which are 
accompanied by a second-order differential equation de- 
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scribing the equation of motion of the perturbed dark en- 
ergy. Even in modified gravity context, in the Hterature, 
we often notice a similar equation replacing G by some 
effective Goff. Without proper (perhaps numerical) veri- 
fication such a simplification is hardly allowed mathemat- 
ically because it corresponds to replacing a second-order 
differential equation by an algebraic coefficient (which 
is zero in the above case); as we have GeS = G in Ein- 
stein's gravity limit, if such an approximation of ignoring 
the DEP is not allowed in Einstein's gravity the same is 
true even in modified gravity context. 

Indeed it is always prudent and correct to include 
the DEP in principle, but more relevant issue would be 
whether we could ignore such accompanied fluctuations 
in practice. In this Letter, by using a simple dynamical 
dark energy model based on a scalar field we will show 
that the answer is negative even in Einstein's gravity; for 
related works, see [1]. 

As a simple dynamical model of dark energy we con- 
sider a minimally coupled scalar field with a double 
exponential potential (we set c = 1 = /i) V{(j)) = 
Fie"^!'^ ^26-^2'^, where is the scalar field. The back- 
ground evolution was investigated previously by Bassett 
et al. in and we consider the background parameters 
of the scalar field suggested in that work: we call it a 
(/)CDM (cold dark matter) model. As a fiducial model we 
take a flat ACDM universe with parameters ri,„ = 0.274 
i^c = 0.2284 and = 0.0456), = 0.726, h = 0.705, 
Us = 0.960, erg = 0.812, To = 2.725 K, Yue = 0.24, 
N„ = 3.04 based on the WMAP 5-year observations [J, 
but without reionization. Evolution of the background 
world models is presented in Fig. [1] For all 0CDM mod- 
els we take Vi = 10~^^ and Ai = 9.43, and from red 
to violet curves, A2 = 1.0, 0.5, 0.0, -0.2, -1.0, -10, 
and —30. For each model, V2 parameter has been deter- 
mined to have the present dark energy density parame- 
ter equal to fi^ = 0.726. The initial dark energy density 
parameter rt^i is determined by the parameter Ai; i.e., 
^(j>i = 3(1 + 11)) /Xf = 0.045 during the radiation domina- 
tion with w = i, see Eq. (4) in [U\ . 

Our dark energy model allows exact scaling during the 
radiation and matter dominated eras (provided by Ai- 
term) and behaves as the dark energy in the present 
epoch (provided by A2-term). Following we consider 
the initial contribution from the dark energy to be close 
to a maximum amount allowed by the bi g ba ng nucle- 
osynthesis (BBN) calculation fi^^ < 0.045 [H]. The pa- 
rameters used in our dark energy model are consistent 
with currently known cosmological constraints from the 
BBN and the high-z SNIa observations, see Fig.[TJ 

In order to calculate the matter and CMB power spec- 
tra, and evolution of the baryon density perturbation we 
solve a system composed of matter (dust and CDM) , ra- 
diation (handled using the Boltzmann equation or tight 
coupling approximation) , together with the cosmological 
constant or the scalar field as the dark energy. Our set 
of equations and the numerical methods are presented 
in [l3|. As the initial conditions for perturbation vari- 
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FIG. 1: Top panels: Evolution of fli and pi as a func- 
tion of scale factor a{t) in the 0CDM universes with scalar 
field potential parameters set by Basset et al. (colored 
curves), where i = T,m,(f> indicates radiation, matter (baryon 
-I- CDM) , and scalar field, respectively. Black curves represent 
those of ACDM model. Middle and bottom panels: Evolution 
of r^^, w^, Hue{z)/Hacdm, and A/i(z) = ^de(z) -Macdm(z) 
for the same set of 0CDM models. In the A^-plot, the grey 
open squares with error bars represent the deviation of SNIa 
data points from the ACDM model considered here. The 
binned SNIa data are based on the Union sample [13] . 



ables we use the scaling solutions derived in |TT|. We 
solved the system in three different gauge conditions: 
the CDM-comoving gauge (CCG), the uniform-expansion 
gauge (UEG), and the uniform-curvature gauge (UCG); 
the CCG, the UEG, and the UCG, respectively, set the 
velocity of the CDM, the perturbed expansion of normal 
frame vector (or the perturbed trace of extrinsic curva- 
ture) , and the perturbed part of intrinsic scalar curvature 
equal to zero as the temporal gauge condition; all pertur- 
bation variables we use are spatially gauge-invariant [TJl ■ 
The CCG is the same as the synchronous gauge without 
the gauge mode. Each of these gauge conditions fixes 
the gauge degrees of freedom completely, thus any vari- 
able in these gauge conditions is equivalent to a unique 
gauge-invariant combination of variables. Value of any 
gauge-invariant variables evaluated in the three gauges 
should coincide exactly. We used this to check the con- 
sistency of the calculation and the numerical accuracy. 

In Fig. [5] we present the matter power spectrum and 
the CMB temperature and polarization anisotropy power 
spectra based on the same parameters used in the back- 
ground world model. The CMB temperature and polar- 
ization anisotropics are naturally gauge-invariant, and for 
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FIG. 2: The matter power spectrum (top-left), and CMB 
TT (top-right), EE (bottom-left), TE (bottom-right) power 
spectra of <^CDM universe with scalar field potential param- 
eters used in Fig. [1] with the same colored code. Predictions 
of ACDM model are shown as black curves. The vertical 
line in the top-left panel indicates the present horizon size 
(10081 Mpc) of the ACDM universe. The small box in 
the top-right panel magnifies the CMB TT powers at low ^'s. 
All calculations are made in three difi^erent gauge conditions 
(CCG, UEG, and UCG), where evolution of perturbation of 
the dark energy scalar field has been properly considered. The 
results in the three gauges coincide exactly. The matter and 
CMB power spectra of the ACDM model have been normal- 
ized with £78 and COBE spectrum, respectively. For com- 
parison, all the (;/)CDM power spectra have been normalized 
with the ACDM ones at small scales, I = 700 for CMB and 
fc = 0.3 /i Mpc"^ for matter ones. For a <t>CQM with A2 = 1.0 
that is most deviated from the ACDM prediction, the ratios 
of its powers to our ACDM predictions are also shown in the 
bottom region of top panels; as an indication of numerical 
accuracy of our code "the CMBFAST-derived power spectra 
[lil divided by our result for ACDM model" is represented as 
black curve. 



the matter power spectrum we present the power spec- 
trum of density perturbation based on the CCG which 
is also a gauge-invariant concept; i.e., density perturba- 
tion in the CCG is the same as a unique gauge-invariant 
combination between the density perturbation and the 
velocity perturbation of the CDM component. Despite 
the variety of outcome in the rcdshift-distance relation 
in the parameters used (sec right-bottom panel in Fig. 
[1]) the matter power spectra of the t/fCDM models are all 
similar with some tilt relative the fiducial ACDM model, 
whereas the differences in the CMB power spectra are 
less distinguished. Figure [2] shows that when we properly 
include the DEP the three gauges give identical results 



X [h-' Mpc] e-180"/i [degree] 

10' 1000 100 100 10 1 0.1 




FIG. 3: The same as Fig. [2] now ignoring the DEP (DEP- 
OFF) in the CCG (blue, dashed), the UEG (green, long 
dashed), and the UCG (brown, dotted curves) for A2 = 1.0. 
Red solid curves represent power spectra with proper DEP 
(DEP-ON). The power spectra ignoring the DEP apparently 
depend on the gauge choice which reflects internal inconsis- 
tency of the system. For matter and CMB TT power spectra, 
recent measurements from SDSS DR7 LRG [3 and WMAP 5- 
year data (including the cosmic variance) have been added 
(grey dots with error bars), and power ratios between cases 
ignoring and considering DEP are also shown for CCG, UEG, 
and UCG, together with fractional errors of observed spectra. 



both for the ACDM and 0CDM cases. 

Now, in Fig. [3] we ignored (set equal to zero by hand) 
the perturbed part of dark energy. Apparently, the re- 
sults depend on the gauge conditions used. As the values 
of gauge-invariant variables depend on the gauge condi- 
tions used in the calculation this alarms inconsistency 
of the system. Such differences are expected because by 
ignoring the DEP the perturbed system of equations be- 
comes inconsistent. The presence of fluctuations in the 
matter and metric simultaneously and inevitably excites 
fluctuations in the dark energy. And it is not allowed to 
turn off the DEP by hand. The issue we would like to 
address, however, is whether we could ignore the DEP 
in practice. Our result in Fig. [3] shows that ignoring the 
DEP easily leads to observationally significant deviations 
in the power spectra which are even gauge dependent. 

In our normalization the matter power spectrum shows 
about -20%/-34%/-|-20% (-10%/-19%/-H8.9%) error 
caused by ignoring the DEP at fc ~ 0.022/iMpc~^ in 
the CCG/UEG/UCG; the values inside parenthesis are 
for Q.tf,i = 0.0225 which is one half of the value used 
in our Figures. The current observation from SDSS 
DR7 LRG (Luminous Red Galaxies) shows 11% (corre- 
lated) error at the same scale, which is already smaller 
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FIG. 4: Evolution of baryon density perturbation (top-left), 
and the normalized perturbation growth factor g = {5b /a) in 
three different scales for A2 = 1.0; due to strong deviations 
we omit UCG cases in top-left panel. As we normalize g to 
unity at present, the effect of DEP appears only in the large 
scale (top-right), and except for the UCG, it has no effects in 
the two small scales (bottom panels). We add 1% error bar 
expected from future X-ray and weak lensing observations [^. 



than the deviations caused by ignoring the DEP in all 
gauges. The CMB temperature power spectrum shows 
about -9.8%/-18%/+.64% (-6.0%/-10%/+.63%) er- 
ror caused by ignoring the DEP at £ = 200; the WMAP 
5-year data in this scale has about 2% (binned) error 
(mostly due to the cosmic variance); as the figure shows 
the small error in the UCG is only due to a coincidence 
in this scale. 

Thus, deviations depend directly on the amount of 17^^ 
in the early scaling era. In our scaling dark energy model 
by reducing the deviations caused by ignoring the 
DEP become proportionally smaller. However, this does 
not imply that our model is an extreme example in the 
effects of DEP. In fact, we can easily introduce models 
theoretically (i.e., not by hand) where is negligible 
during the nucleosynthesis era but becomes significant 
during later radiation and early matter eras and then re- 



duces to the dark energy in recent era so that the rest 
of the cosmological effects are indistinguishable but the 
resulting power spectra in the matter and CMB are sub- 
stantially different in diverse ways, see [l6j . 

This still implies that the substantial deviations in the 
power spectra due to DEP are mainly caused during the 
scaling era. This is partly supported by studying the 
baryon density perturbation growth factor g in the recent 
past which provides another domain where theory meets 
with observation, see Fig. [H In the CCG and the UEG 
cases the DEPs do not cause difference in observation- 
ally relevant small scales. Although the observationally 
distinguishable substantial deviations in the UCG case 
can be regarded as exceptional peculiarity of that gauge 
choice, this still indicates the inconsistency of equations 
without DEP and potential danger of ignoring the DEP 
without proper confirmation. That is, by ignoring the 
DEP the system of equations becomes inconsistent and 
even (gauge-invariant) observable results depend on the 
gauge choice; thus. Fig. 3] shows the particular impor- 
tance of taking proper gauge in the absence of the DEP. 
Notice that in our realistic situation with early radiation 
era the growth factor shows scale dependence. 

In this Letter we investigated the roles of DEP in a dy- 
namical dark energy model based on the scalar field. The 
moral is that when we consider dynamical dark energy it 
is essentially important to take into account the fluctu- 
ating aspects of dark energy properly. When one ignores 
DEP it is important to show that one can do that with- 
out causing observationally significant differences. Our 
model shows an example where it is crucially important 
to include the DEP. Otherwise, the system of equations 
becomes inconsistent, and the consequent results are not 
reliable compared with currently available observations. 
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